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Abstract 

Inclusive cross sections for Higgs boson production in proton-proton 
collisions are calculated in the formalism of unintegrated gluon distri- 
butions (UGDF). Different UGDF from the literature are used. Al- 
though they were constructed in order to describe the HERA deep- 
inelastic scattering Fi data, they lead to surprisingly different results 
for Higgs boson production. We present both two-dimensional in- 
variant cross section as a function of Higgs boson rapidity and trans- 
verse momentum, as well as corresponding projections on rapidity or 
transverse momentum. We quantify the differences between different 
UGD's by applying different cuts on interrelations between transverse 
momentum of Higgs and transverse momenta of both fusing gluons. 
We focus on large rapidity region. The interplay of the gluon-gluon 
fusion and weak-boson fusion in rapidity and transverse momentum is 
discussed. We find that above px ~ 50-100 GeV the weak-gauge-boson 
fusion dominates over gluon-gluon fusion. 

PACS 12.38.Bx,12.38.Cy,13.85.Qk,14.70.Hp,14.80.Bn 



1 Introduction 

Recently unintegrated gluon (parton) distributions became a useful and in- 
tuitive phenomenological language for applications to many high-energy re- 
actions (see e.g. |H 12] and references therein). Mostly the HERA F 2 data 
was used to test or tune different models of UGDF's. However, the structure 
function data are not the best tool to verify UGDF, in particular its depen- 
dence on gluon transverse momentum, because it enters to the j*p total cross 



1 



section in an integrated way. UGDF's have been used recently to describe jet 
correlations [3], correlations in heavy quark photoproduction |[|, total cross 
section for Higgs production [Hj, inclusive spectra of pions in proton-proton 
collisions |H| or even nucleus- nucleus collisions 7j. It is rather obvious that 
differential cross sections seem a much better tool than total or integrated 
cross sections to verify UGDF's. 

Many unintegrated gluon distributions in the literature are ad hoc parametriza- 
tions of different sets of experimental data rather than derived from QCD. 
An example of a more systematic approach, making use of familiar collinear 
distributions can be found in Ref. |Hj. Recently Kwiecihski and collaborators 
pirnHrrTI have shown how to solve so-called CCFM equations by introducing 
unintegrated parton distributions in the space conjugated to the transverse 
momenta . We present results for inclusive Higgs production based on un- 
integrated gluon distributions obtained by solving a set of coupled equations 
jTTj. Recently these parton distributions were tested for inclusive gauge bo- 
son production in proton-antiproton collisions and for charm-anticharm 
correlations in photoproduction jl]. 

While in the gauge boson production one tests mainly quark and anti- 
quark (unintegrated) distributions at scales fi 2 ~ M$y, M|, in the charm- 
quark photoproduction one tests mainly gluon distributions at scales /i 2 ~ 
m 2 c . The nonperturbative aspect of UPDF's can be tested for soft pion pro- 
duction in proton-proton collisions [15] . 

Different ideas based on perturbative and nonperturbative QCD have 
been used in the literature to obtain the unintegrated gluon distributions in 
the small-x region. Since almost all of them were constructed to describe the 
HERA data it is necessary to test these distributions in other high-energy 
processes in order to verify the underlying concepts and/or approximations 
applied. It is the aim of this paper to show predictions of these quite different 
UGDF's for Higgs production at LHC at CERN although we are aware of 
the fact that a real test against future experimental data may be extremely 
difficult. We compare and analyze two-dimensional distributions for inclu- 
sive Higgs production in rapidity and transverse momentum (y,pt) to study 
a potential for such an analysis in the future. We focus not only on midra- 
pidities but also try to understand a potential for studying UGDF's in more 
forward or backward rapidity regions. The results of the gluon-gluon fusion 
are compared with other mechanisms of Higgs boson production such as WW 
fusion for instance. 
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2 Formalism 



There are several mechanisms of the Higgs boson production which have 
been discussed in the literature. Provided that the Higgs mass is larger than 
100 GeV and smaller than 600 GeV the gluon-gluon fusion (see FigJTJ) is the 
dominant mechanism of Higgs boson production at LHC energies [T3|. The 
WW and ZZ fusion is the second important mechanism for the light-Higgs 
scenario. Often the so-called associated Higgs boson production (HW, HZ 
or Htt) is considered as a good candidate for the discovery of the Higgs 
boson. The contribution of the associated production to the inclusive cross 
section is, however, rather small. In the present paper we concentrate on the 
inclusive cross section and in particular on its dependence on rapidity and/or 
transverse momentum of the Higgs boson. 



2.1 Gluon-gluon fusion 

In the leading-order collinear factorization approach the Higgs boson has 
a zero transverse momentum. In the collinear approach the finite trans- 
verse momenta are generated only at next-to-leading order. However, the 
fixed-order approach is not useful for small transverse momenta and rather a 
resummation method must be used jTH]. Recently Kwiecihski, starting from 
the CCFM equation jT3], has proposed a new method of resummation 
based on unintegrated gluon distributions. 

Limiting to small transverse momenta of the Higgs boson, i.e. small 
transverse momenta of the fusing gluons, the on-shell approximation for the 
matrix element seems a good approximation. There are also some technical 
reasons, to be discussed later, to stay with the on-shell approximation. In 
the on-shell approximation for the gg — ► H transition matrix element 1 the 
leading order cross section in the unintegrated gluon distribution formalism 
reads 2 

dof* _ ^gg^H f f / 2 ,,2\ r / „2 ,,2\ X2/-J , -j -> \ d 2 K\ (pKg 



j to — — / Jg/WL: n l: P ) Jg/2\-^2i ™2> P ) u V"-l ">" n 2 ~ ft) 

aya Pt J n n 

(1) 

In the equation above the delta function assures conservation of transverse 
momenta in the gluon-gluon fusion subprocess. The 1/tt factors are due to 
the definition of UGDF's. The momentum fractions should be calculated as 



fg/l(x!, K t , fi ) f g /2(x2,K 2 ,V ) & («1+K2-Pt) 



1 There is a disagreement in the literature ^JEO] how to include the off-shell effects. 

2 Some of UGDF's from the literature depend only on longitudinal momentum fraction 
and transverse momentum of the virtual gluon. To keep formulae below general we allow 
for a scale parameter needed in some distributions. 
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x i,2 — exp(±y), where in comparison to the collinear case M# is replaced 
by the Higgs transverse mass rrit,H- If we neglect transverse momenta and 
perform the following formal substitutions 

fgfifa, >4, V?) ^ X 2 g 2 (x 2 , /J 2 ) <J(«|) . 

we recover the well known leading-order formula 



dyd?p t 



al 9 ^ H x 1 g 1 {xi, ji 2 ) x 2 g 2 (x 2 ,fi 2 ) S 2 (p t ) . (3) 



The off-shell effects could be taken into account by inserting the 7*7* — > 
H off-shell cross section (corresponding to the matrix element squared | Ai"~* H (k,\, k 2 ) 
under the integral in the formula (0) above. This will be discussed in more 
detail in a separate subsubsection. 

There are a few conventions for UGDF in the literature. In the convention 
used throughout the present paper the unintegrated gluon distributions have 
dimension of GeV -2 and fulfil the approximate relation 



f g (x,K (,/x )) dn ^xg coU (x,/j, ) , (4) 







where g co u(x,fi 2 ) is the familiar conventional (integrated) gluon distribution. 
The scale /z 2 in UGDF above is optional. In the effective Lagrangian approx- 
imation and assuming infinitely heavy top quark the cross section parameter 

a o 9 ^ H i s gi ven by 123 

In the following we shall take /1 2 , = M\. Above we have assumed implicitly 
that the fusing gluons are on- mass-shell. In general, the fusing gluons are 
off-mass-shell. This effect was analyzed in detail in the production of Higgs 
associated with two jets [30! • The effect found there is small provided Mh < 
2m t . 

The UGDF's are the main ingredients in evaluating the inclusive cross 
section for Higgs production. Depending on the approach, some UGDF's 
EBl EHj in the literature depend on two variables - longitudinal momentum 
fraction x and transverse momentum k 2 , some in addition depend on a scale 
parameter [HI El EH Ell • In the last case the scale /i 2 is taken here as M\ or 
£M|f, where £ is some factor. 
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The seemingly 4-dimensional integrals in Eq.fJTJ can be written as 2- 
dimensional integrals after a suitable change of variables «i,«2 —> Pt,Qt, 
where pt = K\ + ^2 and qt = k,i — k,2- Then 



where «i = p t /2 + q t /2 and k 2 = p t /2 — q t /2. The integrand of this "reduced" 
2-dimensional integral in g t = ^ - k 2 is generally a smooth function of q t 
and corresponding azimuthal angle <f> qt . 

2.1.1 Unintegrated gluon distributions 

In the present analysis we shall use different unintegrated gluon distributions 
from the literature. We include gluon distributions corresponding to a sim- 
ple saturation model used by Golec-Biernat and Wusthoff to describe the 
HERA deep-inelastic data (GBW), a saturation model of Kharzeev and 
Levin used to describe rapidity distributions of charged particles jZj (KL) 
3 , Balitskij-Fadin-Kuraev-Lipatov (BFKL)-type UGDF and three other 
distributions in the transverse-momentum space j2Zj, [Hj (KMR) and [2E] as 
well as the Kwiecihski UGDF in the b-space. A more detailed description of 
almost all distributions mentioned above can be found e.g. in Ref.JBj. The 
Kwiecihski UGDF as the only one defined in the b-space requires a separate 
discussion. It will be shown below that the formulae for the inclusive cross 
section for Higgs boson production via gluon-gluon fusion can be written in 
the equivalent way also in terms of distributions in the b-space. 

2.1.2 Kwiecinski gluon distribution and the inclusive cross section 

Taking the following representation of the 5 function 



the formula (JBJ) can be written in the equivalent way in terms of gluon dis- 
tributions in the space conjugated to the gluon transverse momentum 4 



3 The normalization of the gluon distributions was fixed in [^] to reproduce the HERA 
data. 

4 The simple form of the formula below would not be possible with off-shell effects, i.e. 
when \A4"^ H | 2 is a function of K\ and K2. 




(6) 
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where 

/■oo 

f g {x,b,n 2 )= dK t K t J {K t b)fg(x,^,fi 2 ). (9) 
Jo 

For most of unintegrated gluon distributions formula © is used. Formula 
((HI) is used when applying the Kwiecihski unintegrated distributions obtained 
as a solution of his equations in the b-space. The b-space approach proposed 
by Kwiecihski is very convenient to introduce the nonperturbative effects like 
intrinsic (nonperturbative) transverse momentum distributions of partons in 
nucleons. It seems reasonable, at least in the first approximation, to include 
the nonperturbative effects in the factorizable way 

f g (x, b, /i 2 ) = j£ GF1£ {x, b, /i 2 ) • F7(6) . (10) 

The form factor responsible for the nonperturbative effects must be normal- 
ized such that 

F» p {b = 0) = 1 . (11) 

Then by construction: 

f g (x,b = 0,n 2 ) = ^g(x,n 2 ) . (12) 
In the following, for simplicity, we use an x-independent form factor 

F7( b ) = «P ('Q (13) 

which is responsible for the nonperturbative effects. The Gaussian form 
factor in b means also a Gaussian initial momentum distribution exp(— k 2 b"l) 
(Fourier transform of a Gaussian function is a Gaussian function). Gaussian 
form factor is often used to correct collinear pQCD calculations for so-called 
intrinsic momenta. Other functional forms in b are also possible. 

The similarities and differences betwen the standard soft gluon resumma- 
tion and the CCFM resummation have been discussed recently in Ref.|17j. 
It has been shown how the soft gluon resummation formulae can be obtained 
as the result of the approximate treatment of the solution of the CCFM 
equation in the so-called b-representation. 

2.1.3 Off-shell matrix element for g*g* — > H 

While in the collinear approach the off-shell matrix elements are needed only 
at higher orders JH] , i n the ^-factorization approach the off-shell matrix ele- 
ments appear in principle already at the leading order. In Ref . the matrix 
element was calculated in the framework of the effective Lagrangian for the 
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Higgs boson coupling to gluons in the infinitely heavy top mass approxima- 
tion. While the on-mass-shell couplings in the full theory and in the effective 
Lagrangian theory are equivalent, this is not expected for matrix element 
with off-shell gluons. In particular the dependence on transverse momenta 
and their relative orientation can be different. The use of the effective vertex 
all over the phase space may not be completely realistic. However, the matrix 
element in the full theory (with finite top mass) was not yet calculated. 

In the present paper we shall only estimate the off-shell effects in the 
effective Lagrangian approximation. Then the cross section for the Higgs 
production can be written as: 

au u It/ 2 2\ e I 2 2\ o 1 ' 2 \t 2 2 2\ j2 

~dy~¥p~ t = l^nf J fgflfa' K i> V ) fg/2{X2, k 2 , fi ) 2 HKi,K2,Pt) d q t 

(14) 

where A is a smooth function of its parameters fOlEHl- In the limit p\ — ► 0, 
k\ — > 0, «2 — > the dimensionless function \(k\,k\,p1) — > 1. Then the 2 
cos 2 (f>j( lt g 2 factor in the formula above constitute the essential difference with 
respect to the on-shell approximation. It modifies the dependence of the 
integrand under the <p q integration. 
In FigElwe show the ratio: 

R = (da /dp t )of /shell / 15 x 

(da /dp t ) on- shell 

as a function of Higgs-boson transverse momentum. In this calculation the 
BFKL unintegrated gluon distributions have been used as an example. The 
result for other distributions is similar. The final effect is rather small (less 
than 10%) in the region of our interest and depends on the way how the 
flux factor for off-shell gluons is defined [T5J 120] 5 - At finite, but not too 
large, Higgs boson transverse momenta the averaging over gluon transverse 
momenta with UGDF's gives < 2cos 2 (f) Kl)R2 >~ 1, and one approximately 
recovers the on-shell result. This is not true for p t w 0, when K\ and are 
strongly anticorrelated and the averaging is not efficient. The off-shell effect 
on the integrated cross section (da/dy) is even smaller. Here we wish to 
concentrate rather on the effect of transverse momenta inherent for UGDF's. 
We shall leave a detailed study of the off-shell effects in effective Lagrangian 
and in the full theory for the future and consequently we shall use the on- 
shell matrix element in the following. This approximation will be also useful 
here when comparing the /^-factorization results with that for the standard 



5 While it is rather strightforward to calculate matrix element for off-shell gluons, there 
is some ambiguity in defining the flux factor for virtual gluons. 
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collinear and soft-gluon resummation approaches. The on-shell approxima- 
tion was used recently in the formalism od doubly unintegrated parton dis- 
tribution for electroweak boson production [T8] . 

2.1.4 Standard soft-gluon resummation 

The formula for inclusive cross section in terms of unintegrated gluon distri- 
butions in the impact parameter space looks very similar to the one in the 
standard soft gluon resummation approach known from the literature. This 
similarity is not random J7j. In the Collins-Soper-Sterman (CSS) approach 
[T6] the resummed cross section for Higgs production reads 



X\ 



%2 



da 



a, 



gg^H 
o 



dyd 2 p t:H (2vr) 



d 2 bJ (p t b)W^(b, Xl ,x 2 ,^) 



h 



2tt 



g 2 (x 2 ,p(b)) + ^Wc vc g 2 (x2,Kb)) + 

Z7T 



2tt 
a s (/j,(b)) 



2tt 



exp 



-(S g (b,^)+S g (b,^)) 
(16)" 



where the exponents in the Sudakov-like form factors read 



S g (b, /i 2 



dfj/ 

(6) P 2 



hi 



A g (a a (fj, 2 )) + B g (a s (fi 2 )) 



(17) 



The coefficient functions C's in Eq. fllfij) can be found in [21]. The coefficient 
A and B in the Sudakov-like form factor can be expanded in the series of a s : 



A„ 



B n 



Ca ^T + 



-2A 



2tt 



+ 



2vr 
2tt 



+ ... 



+ ... 



where /3 = fC A - \N F T R {T R = ±,N F = 5,C A = 3). The CSS formalism 
[T6] leaves open the question of small b. Different prescriptions have been 
proposed to treat this region. The lower limit of the integral in Eq. (jl7j) 

is usually taken /^„ in (6) = (-jf) , where C& = 2 exp(— 7^) « 1.1229. This 
prescription leads to a kink for the Sudakov form factor if Cb/b = \i. To allow 
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a smooth dependence and to quarantee that the lower limit is really lower 
than the upper limit, one could make the following replacement \J? min {ti) = 

(ir) 2 i^r) 2 i 1 + C b / ( b2 V 2 )} 1 ■ To quarantee that the scale of parton 
distribution does not take unphysically small value we shall use the following 
prescription: 

^ 2 (b)=^ mm (b)+ixl, (19) 

where /i 2 , is the starting value for the QCD evolution. In the present paper 
we shall use easy to handle leading order parton distributions from Ref . [2H] • 
W^i p (b, xi, x 2 , /i 2 ) in Eq.ffTtjj) is of nonperturbative origin. Different ef- 
fective parametrizations have been proposed in the literature. Assuming a 
factorizable form of the W^ g p function 

W g N g p (b, xi, a*, fi 2 ) = F g NP (b, x u /i 2 ) " F g NP (b, x 2 , /i 2 ) (20) 

the soft-gluon-resummation formula (fT7)|) and the unintegrated gluon distri- 
bution formula (JSJ) for Higgs production in the b space have identical struc- 
ture if the following formal assignment is made: 

/f^xA^)^/^ • (21) 

If the off-shell matrix element for gg —>■ H is taken in the UGDF ap- 
proach with the Kwiecihski distribution, the structure of the formula in both 
approaches would be different. In this sense the UGDF approach seems more 
general than the b-space resummation method. 

2.2 2^2 processes 

At sufficiently large transverse momenta (p t > Mjj) the Higgs boson produc- 
tion of the type 2 — > 2 should dominate over the 2 — > 1 mechanism discussed 
above. The cross section for fixed-order processes of the type pip 2 — > Hp 3 
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(parton+parton —> Higgs+parton) of the order of a s is well known [31] 
da ,1 



dy H dy p d 2 p, 



ivw,Vp,pt) 



x <^ xi5fi(xi,/i 2 ) x 2 g 2 (x2, V 2 ) \M 



2\ TT~A 12 

gg^Hg\ 



+ 



./i=-3,3 



X 2 g 2 (x 2 , H ) \M qg ^ H q\ 



x 2 q 2 j 2 (x 2 , ii 2 ) 



./a=-3,3 



gq^Hq\ 



^ xigij(xi ;/ u 2 ) x 2 g 2 ,-/(a;2,/i 2 )|A<gq^Hg| 2 

— Q Q ^ 



/=-3,3 



:22i 



The indices / in the formula above number both quarks (/ > 0) and anti- 
quarks (/ < 0). Only three light flavours are included in actual calculations. 
The explicit formulae for \Ai\ 2 can be found in 



2.3 Weak boson fusion 

Up to now we have discussed only the contribution of the dominant LO 
gluon-gluon fusion and NLO 2^2 corrections and completely ignored con- 
tributions of other processes. The second most important mechanism for 
Higgs production is the fusion of off-shell gauge bosons: WW or ZZ. It is 
known that at LHC energy and intermediate mass (100 GeV < Mh < 500 
GeV) Higgs the WW fusion constitutes about 10-15 % of the integrated in- 
clusive cross section. If the weak boson fusion contribution was separated, 
the measurement of the WWH (or ZZH) coupling would be very interesting 
test of the Standard Model. 

Previous studies of the WW mechanism concentrated on the total cross 
section for the Higgs production. In the present paper we are interested in 
differential distributions of Higgs boson rather than in the integrated cross 
section. 

For the gauge boson fusion the partonic subprocess is of the 2-^3 type: 
q(pi) + qfa) q(p3) + qip*) + H(ph)- The corresponding hadronic cross 
section can be written as 

_ VV{ \ ]_Tm 12 d3 Pz d 3 Pi d3 PH 

aa-s 12 [x u x 2) 2§ \Jvi qq ^ qqH \ (27r )3 2 £ 4 (2tt) 3 2^ (23) 

(27r) 4 5 4 (pi +P2-P3-P4- Ph) dx x dx 2 ■ 
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The next-to-leading order corrections to the matrix element of the WW fusion 
are rather small F° r comparison the NLO corrections for gluon-gluon 
fusion are significantly larger. Since we wish to concentrate on relative effect 
of the gluon-gluon and WW fusion contributions in the following we restrict 
to a much simpler leading order (LO) calculation. The LO subprocess matrix 
element was calculated first in Ref . |33j . The spin averaged matrix element 
squared reads 

1 1 F (2p 3 - Pl + M^(2p 4 -p 2 + M^ K ' 

For the WW fusion, limiting to light flavours, the partonic function is 

J- l2 [Xi,X2) — 

(ui(xi,ij,{) + + (u 2 (x 2 ,f4) + d 2 (x2,t4) + s 2 (x 2 ,[il)) + 

(ui(xi,nl) + diix^fil) + si(xi,nl)) (u 2 (x 2 ,nl) + d 2 (x 2 ,f4) + s 2 {x 2 ,n\)) • 

(25) 

We take either (i) fi\ = fi% = M H or (ii) a\ = —ti,/j% = —t 2 , where t\ and t 2 
are virtualities of W bosons. It is convenient to introduce the following new 
variables: 

(26) 

P- = P3 ~ P4 , 

which allow to eliminate the momentum-dependent <5 3 (...) in Eq. (j23j) . Instead 
of integrating over x\ and x 2 we shall integrate over yi = ln(l/xi) and 
y 2 = ln(l/£ 2 )- Then using Eq.(|23|) we can write the inclusive spectrum of 
Higgs as 

da f . n n\ 1 d 3 p- -r-r- TTT 1 1 

— =J dy ld y 2 x lX2 H^x 2 ,^^ - — \M qq ^ — — 
^2^)5 ^ ^ l + E 2 — E 3 — E4 — E H ) ■ 

(27) 

This is effectively a four-dimensional integral which can be easily calculated 
numerically. 

3 Results 

3.1 Gluon-gluon fusion 

Since we wish to concentrate on the potential to verify different UGDF's 
rather than to present the best predictions for LHC experiments we shall 
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consider only one mass of the Higgs boson Mh = 125 GeV as an example. 
This is slightly above the lower limit obtained from the analysis of the LEP 
data jHl]. In addition, this is a mass for which many calculations in the 
literature has been performed recently. Therefore this gives a chance of a 
comparison to the existing results. 

Before we go to the analysis of the two-dimensional spectra of Higgs 
boson produced in proton-proton or proton-antiproton collisions let us show 
the range of gluon longitudinal momentum fraction tested in these proceses. 
In FigElwe present correspondingly x\ and X2 as a function of Higgs boson 
rapidity for a few different values of Higgs transverse momentum. While 
at Tevatron energies (panel (a)) only intermediate and large x's come into 
the game, in collisions at LHC energy (panel (b)) x~ 10 -2 is sampled at 
midrapidity. However, at LHC energy, at rapidities \y\ > 2 one enters the 
region of x > 0.1. Here some of the low-x models of UGDF's may become 
invalid. 

Let us concentrate first on transverse momentum distributions. The dis- 
tribution of Higgs transverse momentum (rapidity integrated) is shown in 
Fig0] In FigElwe present transverse momentum distribution of Higgs boson 
in different bins of rapidity specified in the figure caption. At midrapidity 
(panel (a)) only small x's are sampled. Even here different models from the 
literature give quite different transverse momentum distributions although all 
of them give a reasonable description of the HERA data. The LO soft-gluon 
resummation distribution and the distribution obtained with the Kwiecihski 
unintegrated gluon distribution (bo = 1 GeV -1 ) are very similar with max- 
ima at p t ,H ~ 10 GeV and p t ,H ~ 5 GeV, respectively. The cross section 
for W or Z production is fairly sensitive to the choice of the nonperturba- 
tive form factor [TJ]. In contrast to the gauge boson production the Higgs 
boson production is much less sensitive to the parameter of the Gaussian 
form factor. The results with bo = 0.5 or 2 GeV -1 (not shown here) almost 
coincide with the result for bo = 1 GeV -1 . Therefore in the following in all 
calculations we shall use b = 1 GeV -1 . The BFKL type gluon distributions 
lead to much larger cross section at small Higgs transverse momenta and a 
sizeably larger slope of the ^-distribution. The rapid fall-off of the cross sec- 
tion for the Golec-Biernat-Wusthoff nonperturbative gluon distribution [21] 
(thin dashed curve) demonstrates how important is the perturbative initial 
state radiation in generating larger transverse momenta of the Higgs boson. 
Such effects are not taken into account in Ref . |24j . 

A comment regarding the GBW distribution is here in order. This distri- 
bution was obtained based on a simple dipole parametrization of the HERA 
data inspired by the saturation idea. In addition to the very steep distribu- 
tion in Higgs-boson transverse momentum the scale-independent GBW gluon 
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distribution gives very small total cross section (about half of 1 pb). One 
has to remember, however, that this distribution was constructed in order to 
describe oy* p (Q 2 ) for small photon virtualities. We wish to stress here that 
this distribution is not an universal object. For example in its simplest form 
it fails to describe a!fi p (Q 2 ) for large photon virtualities. The corresponding 
effective gluon distribution defined as 

POO 

xg G Bw(x) = / ^ BW (x)dK 2 (28) 

J 

resembles the standard collinear distribution xgnGLAp{x, /i 2 ) for small fac- 
torization scale fi 2 ~ 1 GeV 2 . The latter, when substituted into the standard 
leading-order formula, also leads to a small total cross section of the order 
of 1 pb. A reasonable total cross section is obtained provided /i 2 ~ Mfj. 
Although gDGLAp(x, fi 2 = lGeV 2 ) leads to a reasonable description of F 2 
at Q 2 ~ 1 GeV 2 it cannot be directly (without DGLAP evolution) used for 
Higgs production. In this context the scale-independent GBW distribution 
should be understood as an initial condition for QCD evolution rather than 
an universal object to be used in different high-energy processes. 

Different UGDF's constructed in order to describe the total cross section 
for 7*p process give quite different predictions for the Higgs production. This, 
as will be discussed below, is not completely surprising. While the YP total 
cross section is sensitive to the integral J dk 2 f g (x, fc 2 (,/i 2 )), the Higgs boson 
transverse momentum distribution samples details of the unintegrated gluon 
distributions. As discussed in the previous section, in the case of Higgs boson 
production the inclusive cross section is a convolution of two UGDF's. In 
general, exclusive reactions are much better place for testing UGDF's. In 
this sense the standard procedure to constrain UGDF's through describing 
the F 2 HERA data does not seem very effective. 

In order to understand the situation somewhat better in FigElwe show the 
corresponding average values of sampled transverse momenta as the function 
of Higgs transverse momentum. At y=0, by symmetry requirement, < K\ > 
and < k 2 > are identical. It is not the case for y=3. Rather similar results are 
obtained with different UGDF's. At very small Higgs transverse momenta 
one tests k's ~ 1 GeV. At large Higgs transverse momenta and y=0 we get 
< K\ > + < n 2 >~ Pt/2. Of course by symmetry < Kx/ 2 > (—y) =< k 2 /i > 

(y)- 

Some examples of inclusive rapidity distributions are shown in Fig|7| 
Even here the differences between different UGDF's are clearly visible. Above 
\y\ > 3 only resummation distribution (thin solid line), Kwiecihski distribu- 
tion (thick solid line) and Kimber-Martin-Ryskin (dotted line) are applicable 
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by construction. The other distributions were obtained by extending the gen- 
erally small-x gluon distributions above x > 0.1 by multiplying the original 
formulae by (1 — x) n . The power n=5-7 was found recently in the production 
of cc pairs in photon-proton scattering at low energies jl] . In the present pa- 
per we used n=7. Small differences may be expected only at the very edges 
of the phase space, i.e. in the region we are not interested here. 

For completeness in Fig|H]we show the average values of gluon transverse 
momenta < Ki > and < K2 > as a function of Higgs-boson rapidity. At 
midrapidity by symmetry < K\ >~< K2 >. The average values strongly 
depend on UGDF used and on the region of rapidity. The asymmetry of 
average values of the transverse momenta at forward/backward regions are 
closely related to Fig|3]due to correlation of transverse momenta with longitu- 
dinal momentum fractions. Generally, the smaller x\ {x<i) the larger < K\ > 
(< K2 >)• The details depend, however, on the specific version of the gluon 
dynamics. The extremely small average transverse momenta for the GBW 
UGDF can be understood in the light of the discussion above. 

How important is the choice of the factorization scale in our calculations 
with Kwiecihski UGDF? In FigEl we show results obtained with quite dif- 
ferent choices of factorization scale and with 60 = lGeV -1 . There is very 
little effect if the factorization scale is increased from our cannonical value 
fi 2 = Mjj. There is, however, a sizeable effect if the factorization scale is 
decreased drastically, especially at p t > 50 GeV. The factorization scale de- 
pendence in our case seems somewhat weaker than in a recent work [20J. 

It is particularly interesting to compare the results obtained with the 
Kwiecihski unintegrated gluon distributions with those obtained from the 
standard soft-gluon resummation method. In FigEJJwe show two-dimensional 
distributions in (y,pt)- The distribution obtained with the Kwiecihski UGDF 
decrease less rapidly with Higgs transverse momentum than the distribution 
obtained in the standard soft-gluon resummation. This is partially due to 
the different choice of the factorization scale in both methods. 

A more detailed comparison is made in FigE] where we have selected 
two rapidities y — and y — 3. While at y = the result obtained with 
Kwiecihski distributions and that obtained within the standard resumma- 
tion method almost coincide, they become quite different at y = 3 and p t > 
60 GeV. However, the applicability of both methods at large transverse mo- 
menta is not obvious. The result for y=3 at large transverse momenta ob- 
tained with Kwiecihski UGDF seems more trustworthy than that obtained 
within standard resummation method. The problems of the standard re- 
summation method at large transverse momenta may be caused by some 
somewhat arbitrary prescriptions used as discussed in section 2. 

In order to better emphasize the differences between different UGDF's we 
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separate the contributions to integral in Eq.© from four different disjoint 
and complementary kinematic regions: 

• (I) Ki < p t and k 2 < p t , 

• (II) Ki < p t and k 2 > p t , 

• (III) Ki > p t and k 2 < p t , 

• (IV) Ki > p t and k 2 > p t , 

where k± and k 2 are transverse momenta of the last gluons in the ladders and 
Pt is transverse momentum of the produced Higgs boson. In FigJT2]we present 
the decomposition of the Higgs cross section ^ into those four regions as a 
function of Higgs transverse momentum. Here we limit to midrapidities (-1 
< y < 1) only. In the case of Kwiecihski UGDF, first the Fourier transform 
from the b-space to K t -space was calculated and the results were stored on the 
grid in x and fcf . The grid was used then for interpolation when using formula 
© with the extra conditions on transverse momenta specified above. It is 
interesting to note that at larger transverse momenta the contributions from 
regions II, III and IV are completely negligible. The other contributions are 
important only at very low transverse momenta. For completeness in FiglTTfl 
we show similar decomposition as a function of rapidity. In this calculation 
we have limited Higgs transverse momenta to p t < 40 GeV. The contribution 
of region I dominates at midrapidities. The contribution of asymmetric (in 
K\ and k 2 ) regions II and III becomes important at very forward or very 
backward Higgs production. The contribution of the region IV is negligible 
almost everywhere. The proportions of contributions corresponding to the 
four specified above regions differ significantly for different UGDF's. 

3.2 Estimates of higher-order effects 

In the present paper we have limited to leading-order approach only. This 
was dictated by the fact that until now only leading-order approach was 
used to describe the HERA data in terms of UGDF's. Furthermore the 
consistent next-to-leading order analysis is rather complicated. We leave the 
next-to-leading order analysis for a future study. Instead, we wish to visualize 
(estimate) the NLO corrections in a similar soft-gluon-resummation approach 
where the relevant formalism was worked out in detail |21j . In Fig|T4l we 
present the relevant soft-gluon-resummation K-factor as a function of Higgs 
transverse momentum for -0.5 < y < 0.5 (panel a) and 2.5 < y < 3.5 (panel 
b). In this calculation the Gaussian form factor (see Eq. (fT3| ) with 6 = 1 
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GeV -1 was used. The dashed line includes only gluonic NLO corrections 
to partonic function, the dotted line exclusively the quark NLO corrections 
and the solid line includes all NLO effects altogether. At midrapidities the 
gluonic effects are absolutely dominant and enhance the leading-order cross 
section by about 80 %. The quark corrections are at the level of 1 % and can 
be numerically neglected. They become sizable (of the order of 10 %) at very 
forward and very backward rapidities. In Fig^Jwe present a corresponding 
K-factor (K = NLO/LO) as a function of Higgs rapidity. In this calculations 
Pt,H < 80 GeV was assumed. Summarizing, to a good approximation the 
NLO soft-gluon resummation corrections result in multiplying the LO cross 
section for Higgs production by a factor of about 1.8. 

Up to now we have concentrated at relatively small Higgs transverse mo- 
menta. At high transverse momenta the standard 2 — > 2 mechanisms take 
over. The cross section for perturbative 2^2 processes with Higgs in the 
final state is shown in FiglTHl The gg — > Hg dominates at small transverse 
momenta and midrapidities. The qg — > Hq and gq — > Hq become compa- 
rable to the first contribution at large transverse momenta and forward and 
backward regions, respectively. The contribution of qq — > Hg is negligible 
all over the interesting part of the phase space. 

3.3 Weak boson fusion versus gluon-gluon fusion 

The weak boson fusion is known to be another important ingredient in total 
(integrated) cross section for Higgs boson production [Hj. It is interesting 
what is interrelation between the two dominant contributions in rapidity and 
transverse momentum of the Higgs boson. In FigElwe present such a two- 
dimensional spectrum. This spectrum is very different from those for the 
gluon-gluon fusion. In particular, the maximum of the cross section at p t 
slightly larger than 50 GeV is visible. It is interesting if the contribution of 
weak-boson fusion can exceed in some corner of the phase space the gluon- 
gluon contribution. In order to quantify the effect in FigEl we present d ^ 
as a function of Higgs transverse momentum for y = and y = ± 3. The 
contribution of the WW fusion falls off much faster for y — ± 3 than for 
y = 0. The results are almost independent of the choice of the factoriza- 
tion scale. These results almost coincide (compare thick solid (prescription 
(i) above) and overlaping thick dashed (prescription (ii) above) lines). For 
comparison we present a few examples of the gluon-gluon fusion with BFKL 
(dash-dotted) and Kwiecihski (thin solid) unintegrated gluon distributions 
and standard resummation method with Gaussian form factor and bo = 1 
GeV -1 (dashed). We find that depending slightly on UGDF and rapidity, 
above pt ~ 50-100 GeV the WW fusion mechanism dominates over gluon- 



16 



gluon fusion mechanism. However, the 2^2 processes (dotted line) are 
large, especially at midrapidities. Only at large rapidities and p t > 150 GeV 
the WW fusion seems to dominate over the other processes. However, there 
the cross section is rather small. Whether this opens a possibility to study 
the WWH and similarly ZZH couplings requires further studies. 

4 Conclusions 

In the present paper we have presented predictions for the inclusive cross 
section for the Higgs boson production at the LHC energy W = 14 TeV, 
obtained with the help of different models of unintegrated gluon distribu- 
tions used recently in the literature. Almost all the UGDF's discussed here 
were obtained based on the analysis of low-x HERA data for virtual photon 
- proton total cross sections. Although they are almost equivalent in the 
description of the HERA data, quite different results have been obtained for 
Higgs production. While the structure function data is sensitive to rather 
low transverse momenta of the gluon, in the Higgs production, in principle, 
one could sample not yet explored region of large transverse momenta (at 
large scales). One should remember that the Higgs production even at large 
LHC energy is not completely a small-x phenomenon. The analysis of very 
forward or very backward Higgs boson production, in principle, could open 
a possibility to study UGDF's in a completely unexplored region of large x. 
This task, however, is by no means easy as in this region of phase space both 
"small-x" and "large-x" physics is intertangled. 

We have shown that for all UGDF's discussed here the inclusive cross 
section is dominated by the configurations with transverse momenta of the 
gluons smaller than the transverse momentum of the Higgs ( k± < pt and 
1^2 < Pt )■ 

Finally we wish to emphasize that all the potential to study UGDF in the 
Higgs production discussed here is at present only conditional as it implicitly 
assumes the existence, discovery and good identification of the Higgs boson 
in the future experiments at LHC. We do not need to mention that all this 
would not be possible if the Higgsless scenarios (see e.g. [35] and references 
therein) turned out to be true. Even if the Higgs boson is discovered at LHC 
the statistics may not be sufficient for precise tests of UGDF. 
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Figure 1: Dominant leading-order diagram for inclusive Higgs production for 
Pt < M H . 
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Figure 2: The ratio R as a function of Higgs transverse momentum for the 
full range of Higgs rapidity. The thick dashed line corresponds to neglecting 
the function A, i.e. assuming A = 1. The solid line was obtained with the 
Hautmann prescription of the flux factor ^1 while the dash-dotted line is 
based on the formula from Ref. |2Uj. In this calculation the BFKL UGDF 
was used as an example. 



22 




Figure 3: X\ and X2 as a function of Higgs rapidity y for p t = 0. In panel 
(a) for Tevatron energy \/s = 2 TeV and in panel (b) for LHC energy y/s = 
14 TeV. 
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Figure 4: Transverse momentum distribution of Higgs boson at LHC energy 
W = 14 TeV for different UGDF's from the literature: solid - Kwiecinski, 
thick dashed - KL, thin dashed - GBW, dash-dotted - BFKL. 
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Figure 5: Transverse momentum distribution of Higgs boson at LHC energy 
W = 14 TeV and y = (left panel) and y = ± 3 (right panel) for different 
UGDF's from the literature. The notation here is the same as in the previous 
figure. In addition to the previous figure we present results for KMR ( dotted), 
Kutak-Stasto (grey dash-dotted) and Bluemlein (thin solid) UGDF. 
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Figure 6: Average gluon transverse momentum as a function of Higgs trans- 
verse momentum at LHC energy W = 14 TeV for different UGDF's from the 
literature. 
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Figure 7: Rapidity distribution of Higgs at LHC energy W — 14 TeV for 
different UGDF's from the literature. 
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Figure 8: Average gluon transverse momentum < k\ > or < K2 > as a 
function of Higgs rapidity for different UGDF's from the literature. In this 
calculation p t < 40 GeV. 
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Figure 9: ^4 dependence of the result on the choice of the factorization scale 
for the Kwiecinski UGDF for y = (left panel) and y = 3 (right panel). In 
this calculation a Gaussian form factor with b = lGeV^ 1 was used. 
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Figure 10: A comparison of two-dimensional distributions of Higgs boson for 
(a) Kwiecinski UGDF, (b) LO b-space resummation. 
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Figure 11: A comparison of Higgs transverse momentum distribution calcu- 
lated with Kwiecinski UGDF (red, solid) and LO b-space resummation (black, 
dashed) for different rapidities: y = (left panel) and y— 3 (right panel). 
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Figure 12: Decomposition of the transverse momentum distribution of Higgs 
boson at LHC energy W = 14 TeV and -1 < y < 1 into four regions specified 
in the text. The thick solid line is a sum of all 4 contributions, thin solid - 
region I, dashed - region II+III) and dotted - region IV. 
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Figure 13: Decomposition of the rapidity distribution of Higgs boson at LHC 
energy W — 14 TeV and into four regions specified in the text. In this calcu- 
lation pt < 40 GeV. The thick solid line is a sum of all 4 contributions, thin 
solid - region I, dashed - region II or III) and dotted - region IV. 
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Figure 14: iVXO K-f actor in the soft- gluon-resummation formalism as a func- 
tion of Higgs transverse momentum for two different bins of rapidity. The 
dashed line corresponds to separated gluonic contributions, while the dotted 
line to separated quarkish contributions. The solid line is for both effects 
included simultaneously. 
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Figure 15: NLO K-f actor in the soft- gluon-resummation formalism as a func- 
tion of Higgs rapidity. In this calculation p t < 80 GeV. The meaning of the 
curves here is the same as in the previous figure. 
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Figure 16: Contributions of different subprocesses of the 2-^2 type for W 
= 14 TeV, respectively for (a) gg, (b) qq, (c) qg and (d) gq. 
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Figure 18: Transverse momentum distribution of Higgs boson at LHC energy 
W = 14 TeV and y = (left panel) and y = ± 3 (right panel) produced in 
WW fusion (overlaping thick solid and thick dashed for the two different pre- 
scriptions specified in the formalism section), compared to the corresponding 
contribution of gluon-gluon fusion: BFKL (dash- dotted), Kwiecinski (thin 
solid) UGDF, LO soft-gluon resummation (dashed) and the perturbative 2 — > 
2 collinear contribution (thin dotted). 
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